Rejuvenation and Memory in Migdal-Kadanoff Spin Glasses 
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Abstract. We study aging phenomena of Migdal-Kadanoff spin glasses in order to clarify relevancy of temperature chaos 
to rejuvenation and memory. By exploiting renormalization, we do efficient dynamical simulations in very wide time/length 
scales including the chaos length. As a consequence, we find that temperature chaos and temperature dependence of speed of 
equilibration cause two significantly different effects against temperature variations, i.e., rejuvenation for positive temperature 
variation and memory for negative temperature variation, as are observed experimentally in spin glasses. 



It is well known that in spin glasses, dynamical ef- 
fects strongly depend on the history of the system after 
quench from above the transition temperature T c . These 
phenomena are called aging and have been studied using 
various experimental protocols [1, 2, 3]. Measurement of 
ac-susceptibility during T-cycle [4], which is employed 
in this work, is one of them. This experiment consists of 
the following three stages. In the first stage, the system is 
quenched from above T c and it is kept at a temperature T 
(< T c ) during a time t\. In the second stage the temper- 
ature is temporarily changed to T ±AT (< T c ) during a 
time ?2, and then it is set back to T in the third stage. Ac- 
susceptibility x is measured during all the three stages. In 
the case — AT, % in the third stage resumes its relaxation 
from the value at the end of the period t \ as if the system 
remembers how far the relaxation at T had proceeded be- 
fore the perturbation (memory effect). On the other hand, 
the system is rejuvenated by positive T-cycle, and we 
observe strong relaxation of % in the third stage. 

From a theoretical point of view, "temperature 
chaos" [5], decorrelation of the equilibrium states at 
two temperatures beyond the so-called chaos length 
£(T,T'), has been one of the most conceivable causes of 
rejuvenation. However, temperature chaos seems to be 
incompatible with memory if one naively thinks. There- 
fore, relation among temperature chaos, rejuvenation 
and memory has been vigorously studied [6, 3] as a key 
to understand aging phenomena. In this work, we also 
address this issue by studying Migdal-Kadanoff (MK) 
spin glasses. There are mainly two merits in working 
on MK spin glasses. First the existence of tempera- 
ture chaos is shown in this model, and whose exact 
renormalization [7] allows one to measure the chaos 
length £(T,T'). Second we can do efficient dynamical 
simulations at very long time scales by exploiting renor- 
malization [8]. (See also [9] as a similar approach.) In 
the previous work [8], we have investigated dynamics 
of MK spin glasses by utilizing these advantages, and 
found that temperature chaos causes rejuvenation, but 



it also destroys most of the memory if the length scale 
equilibrated during the second stage is much larger than 
£(T,T'). The main purpose of this work is to show that 
memory is preserved if we take temperature dependence 
of speed of equilibration into account. 

The outline of the paper is as follows. First, we intro- 
duce MK spin glasses. Second, we briefly recall temper- 
ature chaos in MK spin glasses. Third, we explain how 
renormalized dynamics can be used to probe rejuvena- 
tion and memory on very wide time/length scales. Fi- 
nally, we show results of the both positive and negative 
T-cycle simulations. The last section is devoted for dis- 
cussion and conclusions. 

The model — . We consider MK lattices following 
the standard real space renormalization group approxi- 
mation [7] to the Edwards Anderson (EA) model [10]. 
The recursive construction of such hierarchical lattices 
is described in Fig. 1; edges are replaced by 2b edges 
so the "length" of the lattice is multiplied by 2. We call 
generation "level" the order of the recursion and G the 
total number of these. Then the lattice length L is 2 G and 
the number of bonds is (2b) G (which is also roughly the 
number of sites); one can thus identify 1 + ln/?/ln2 with 
the dimension of space on such a lattice. When all the 
edges are constructed, each is assigned a random cou- 
pling jfj. The superscript "0" implies that are bare 
couplings. Similarly, on each site i we put an Ising spin 
Si = ±1. The Hamiltonian is 

Hj({Si}) = - £ jfjSiSj, (1) 

<«/> 

where the sum is over all the nearest neighbor spins of 
the lattice. All of the work presented here will be for 
three dimensions (b — 4) with couplings jfj taken from 
a Gaussian of mean and width 1. The model then 
undergoes a spin glass transition at T c « 0.896 [11]. 

Temperature chaos in MK spin glasses — . A great 
advantage of MK lattices is that it allows us to do renor- 
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FIGURE 1. Construction of a MK lattice (b = 2). 

TABLE 1. Size dependence of C(L, T = 0.7, T 1 = 
0.65). 



L 


2 5 


2 7 


2 9 


2 11 


2 13 


2 15 


C 


1.00 


0.99 


0.92 


0.62 


0.10 


0.00 



malization procedure exactly [7]. Now let us denote a set 
of spins of level n as {S„}. By tracing over spins of lower 
levels, we explicitly obtain 



P k ({So},---,{S G -k})= 



Tr { s G _ t+1 }...{s G }exp(-^°) 

Tr {So}-{s G } ex P(-^°) 
exp[-^/({S },---,{S G -*})],(2) 
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(3) 



where the sum £<(/>„ i s over all the nearest neighbor 
spins of a MK lattice with generation n. Effective cou- 
plings Jfj are Jy/T for k = 0. Otherwise, they are calcu- 
lated by the recursion formula 



jf? 1 = arctanh[tanh^tanhjf 



i=i 



(4) 



In this equation, Jf t and Jfj lie on the l-th path connecting 
i and j. 

Now let us turn to temperature chaos in MK spin 
glasses. If we start from the same set of bare couplings, 
{Jfj{T)} and {Jfj(T + AT)} are completely correlated 
for any T and AT. Does the correlation still survive after 
renormalization is repeated again and again? This ques- 
tion was first addressed by Banavar and Bray [12], and 
they have found that for arbitrarily small Ar {Jfj(T)} 
and {Jjj(T + AT)} become completely decorrelated 
when k is large enough, indicating that spin polarizations 
at two temperatures are different at all. In table 1, we 
show the linear correlation coefficient 



C(L,T,T') 



J k (T)J k (T>) 



(5) 



for T = 0.7 and T = 0.65. In this definition, L = 2 k , ~ 
is the disorder average, a is the standard deviation of J k , 
and we have used the fact that J 1 = 0. We see C(L, T, T 1 ) 
rapidly drops from 1 to around L w 2 U . We hereafter 
utilize C(L, T, T') as an indicator of temperature chaos, 
and define the chaos length £(T,T') as the value of L 
where C = 1/e. (In this case, £(T, T') « 2 12 .) 



Exploiting renormalization for dynamics — . Since 
our purpose is to examine relevancy of temperature chaos 
on dynamics, we should do dynamical simulations at 
very long time so that the length scale equilibrated during 
simulation is comparable with the chaos length. How- 
ever, this condition is hardly satisfied by usual Monte- 
Carlo simulations because frustration and randomness 
inherent in spin glasses make their dynamics extremely 
slow. In fact, measurements of L eq (f), the equilibrated 
length during f, in 3d-EA spin glass model have shown 
that L eq (f) for t = 10 6 Monte Carlo Sweeps (MCS) is 
less than 10 at any temperatures below T c [13], while it 
is almost impossible to go beyond the time scale by the 
present computer resource. This length scale seems to 
be hopelessly shorter than the chaos length. (Recall that 
£(T,T') « 2 12 in the previous case though the tempera- 
ture difference is not so small.) 

In order to overcome the difficulty, we exploit renor- 
malization for dynamics [8]. The basic idea is as follows. 
Suppose we focus on a time window f m ; n < t < f max . Be- 
tween t = and t = f m in the system has had time to equi- 
librate up to the length scale /(f m ; n ); essentially all out 
of equilibrium physics comes from larger length scales. 
On MK lattices, this means that the spins whose gen- 
eration "level" is larger than G m ; n (with 2 NRG = /(f m in) 
and NRG = G G m ; n ) are in local equilibrium; the other 
spins have dynamics that is well described by the effec- 
tive Hamiltonian ^ NRG ({5 }, • • • , {5 Gmin }). In practice, 
we adopt the following procedure to implement this idea 
with taking temperature dependence of speed of equili- 
bration into account. 

1. Calculate the effective couplings at 7h (the higher 
temperature used in T-cycle protocol) and those at 
7l (the lower one). We first generate a large number 
of bare couplings from a Gaussian of mean and 
width 1 . Then, we do renormalization procedure to 
produce an ensemble of effective couplings. This 
process is iterated NRG times. The final effective 
couplings are then randomly assigned to the edges 
of a MK lattice of size 2 Gmin . 

2. The direction of each spin at t = is chosen ran- 
domly with equal probability, corresponding to a 
quench from an infinitely high temperature at an in- 
finite rate. 

3. At 7l, we simply do standard Monte Carlo by using 
^j NRG ({S Q }, • • • , {S Gm J) prepared at step 1. 

4. At 7h, dynamics is further accelerated by 
the following procedure. We first calculate 

JZf RG \{So},---AS G > m J), where NRG' = 
G - G' min and NRG' > NRG. Then we do Monte 
Carlo by using J$?j NRG to update the spins whose 
level is smaller (or equal to) G' min . After each MCS, 
the lower spins {S^} (G' min < k < G m ; n ) are locally 
equilibrated with fixing the spins of smaller levels. 

Note that one MCS at T L (T H ) on the renormalized lattice 
corresponds to a huge number of sweeps on the non- 
renormalized lattice, in fact to the number needed to 



equilibrate on the length scale 2 NRG (2 NRG ). 

Results — . We use the standard temperature cy- 
cling protocol and measure a quantity similar to the ac- 
susceptibility defined as [14] 



l-Q(t- 



2k 



(6) 



where Q(t,t') = L(^'5 ; (f )5 ; (f'))/(L^') and W; is the 
weight of site i. Our results concentrate on the choice 
that Wi is proportional to the connectivity of site i. How- 
ever, our additional simulations of the case Wi = 1 have 
shown that the both choices lead to the same results qual- 
itatively. Every MCS updates all the spins once. An al- 
ternative choice is to sweep the bonds, updating their 
end spins as in [15]; we have checked that the results 
are hardly affected by the method used. In the both pos- 
itive and negative T-cycle simulations, we use Tn = 0.7 
and 7L = 0.65. Note that these two temperatures are the 
same as those in Table 1. The period 2n/(0 of ac-field 
is 16 MCS. All the simulations are done on MK lattices 
with five generations (G m m = 5) using < NRG < 15. 
Since we calculate renormalized couplings at 7r and Tn 
from the same set of bare couplings, they are highly cor- 
related when NRG is small. However, their correlation 
vanishes for large NRG due to temperature chaos. The 
difference NRG' - NRG is 1 for all the simulations. We 
hereafter denote % with a negative (positive) T-cycle as 
^Ncycie CtPcycie) and the isothermal % at T as Xisa(T). 

In Fig. 2, we show three typical behaviors observed 
in positive T-cycle simulations. In the main frame, we 
omit t2 part of data and connect t\ and tj, parts to com- 
pare with the isothermal data drawn by line. For small 
NRG, ^pcycie is remarkably below Xi S o(7L) in the third 
stage, as illustrated in the main frame of Fig. 2(a). This 
means that equilibration at 7l is sharply accelerated in 
the second stage because renormalized couplings at 7l 
and Tii are strongly correlated and equilibration is ac- 
celerated at Tn. This trend (acceleration of equilibration) 
arises until NRG ~ 8. Fig. 2(b) shows that ^p cyc ie be- 
gins to have a strong curvature in the third stage as a 
sign of rejuvenation, while renormalized couplings at 7r 
and Tn are still highly correlated. (As shown in Table 1, 
C = 0.92 for NRG = 9.) However, £p cyc ie is still below 
Ziso(TL) at later times of the third stage. Finally, Fig. 2(c) 
is the result for NRG =13. Renormalized couplings are 
very decorrelated now (C = 0.10). We see strong rejuve- 
nation in the third stage, as is found in experiments. We 
have checked that rejuvenation is perfect in the sense that 
ZPcycie in the first stage and that in the third stage com- 
pletely overlap. We have also found that %p C yc\s in the 
second stage and %i so {Tn) collapse into a single curve 
for any NRGs, meaning that aging at lower temperature 
is not helpful in equilibration at higher temperature. 

For negative T-cycle case, we only show two extreme 
cases in Fig. 3 , i.e., a highly correlated case and a highly 
decorrelated one. Now a surprising fact is that perfect 
memory is observed in the both cases. Especially, in the 
decorrelated case (NRG = 13), memory appears in the 
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FIGURE 2. ^(o),f) with a positive T-cycle (crosses) for 
NRG = 5, 9 and 13 (inset). Temperature is temporarily in- 
creased from 7l = 0.65 to Tn = 0.7 during the period ?2- In the 
main frame, t\ and t^ parts are connected after omitting ?2 part 
to compare with the isothermal data at 7L (line). The average is 
from 6 x 10 3 samples. 



third stage though rejuvenation is observed for positive 
T-cycle (Fig. 2(c)). This result is very contrast with [8] 
which has shown that negative T-cycle leads to strong 
(but not perfect) rejuvenation when NRG = NRG'. These 
findings tell us that temperature dependence of speed 
of equilibration is crucial for memory effect. Lastly, we 
have compared ^Ncycie in the second stage with ^p cyc i e 
in the third stage, and found that they perfectly overlap 
for any NRGs. Since their difference is whether previous 
aging at 7L exists or not, this result suggests that equi- 
libration at higher temperature makes previous aging at 
lower temperature insignificant. 

Discussion and conclusions — . In this work, we 
have found that positive T-cycle and negative one cause 
quite different effect, i.e., rejuvenation in the former and 
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Furthermore, the same behavior is also observed in spin 
glasses if the equilibrated length scale during the second 
stage is much smaller than £(T, T'). (See Fig. 14 of [3].) 
These results may suggest that temperature chaos is ab- 
sent, or the equilibrated length scale in experimental time 
scale is much smaller than £(T, T') for these systems. 
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f3 parts are connected after omitting ti part to compare with 
the isothermal data at 7^ (line). The average is from 6 x 10 3 
samples. 
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